WORLD INTELLECTUAL PROPERTY ORGANIZATION 
International Bureau 




INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(51) International Patent Classification 4 : 
G01R 35/00 


Al 


(11) International Publication Number: WO 89/ 04001 
(43) International Publication Date: 5 May 1989 (05.05.89) 


(21) International Application Number: PCT/US88/03772 

(22) International Filing Date : 25 October 1 988 (25. 10.88) 

(31) Priority Application Number: 1 13,1 10 

(32) Priority Date: 27 October 1987 (27, 10.87) 

(33) Priority Country: US 

(71) Applicant: MARTIN MARIETTA CORPORATION 

[US/US]; 6801 Rockledge Drive, Bethesda, MD 20817 
(US). 

(72) Inventors: SEQUEIRA, Hermann, Brian ; 227 Mallow 

Hill Road, Baltimore, MD 21229 (US). TRIPPE, Mi- 
chael, Wayne ; 1808 Sutton Avenue, Baltimore, MD 
21227 (US). JAKHETE, Rajendra, Sitaram ; 17 Half 
Penney Lane, Catonsville, MD 21228 (US). 


(74) Agents: CHIN, Gay; Martin Marietta Corporation, 
6801 Rockledge Drive, Bethesda, MD 20817 (US) et 
al. 

(81) Designated States: AT (European patent), BE (Euro- 
pean patent), CH (European patent), DE (European 
patent), FR (European patent), GB (European pa- 
tent), IT (European patent), JP, LU (European pa- 
tent), NL (European patent), SE (European patent). 

Published 

With international search report. 
Before the expiration of the time limit for amending the 
claims and to be republished in the event of the receipt 
of amendments. 



(54) Title: CALIBRATING METHOD FOR APPARATUS EVALUATING MICROWAVE/MILLIMETER WAVE 
CIRCUITS 




A calibration procedure for vector network analyzers utilizing a known standard, an unknown standard and a pair 
of offsets bearing a known length ratio (2:1) between them whereby measurements taken on the standards and their combi- 
nations with the offsets result in self-verifying redundant equations which furnish: (1) error terms according to well-known 
flow graph models, (2) the reflection of the unknown standard, (3) the transmission factor of the waveguiding medium 
from which the offsets are realized regardless of the waveguide's loss or dispersion, and (4) a figure of merit (corruption 
coefficient) for the quality of the acquired raw data without the necessity for computing the error terms, connecting verifi- 
cation standards or otherwise completing the calibration process. 



FOR THE PURPOSES OF INFORMATION ONLY 



Codes^edtoidentifyStatespartytothe 



cations under the PCX 



AT Austria 

AU Australia 

BB Barbados 

BE Belgium 

BG Bulgaria 

BJ Benin 

BR Brazil 

CF Central African Republic 

CG Congo 

CH Switzerland 

CM Cameroon 

DE Germany, Federal Republic of 

DK Denmark 

H Finland 



FR 


France 


GA 


Gabon 


GB 


United Kingdom 


HU 


Hungary 


IT 


Italy 


JP 


Japan 


KP 


Democratic People's Republic 




of Korea 


KR 


Republic of Korea 


LI 


Liechtenstein 


LK 


Sri Lanka 


LU 


Luxembourg 


MC 


Monaco 


MG 


Madagascar 



ML Mali 

MR Mauritania 

MW Malawi 

NL Netherlands 

NO Norway 

RO Romania 

SD Sudan 

SE Sweden 

SN Senegal 

SU Soviet Union 

TD Chad 

TG Togo 

US United States of America 



WO 89/04001 



1 



PCT/US88/03772 



CALIBRATION METHOD FOR APPARATUS EVALUATING 
MICROWAVE /MILLIMETER WAVE CIRCUITS 

Background of the Invention 

This invention relates to a calibration method and 
more particularly to a method for calibrating apparatus, 
such as. a vector network analyzer, utilized for making 
5 evaluations of circuits operating at microwave and 
millimeter wave frequencies. 

Most calibration methods used for electronic 
measurement apparatus require the use of a set of extremely 
high quality standards whose characteristics are known in 

10 advance. Additionally, a well known standard must be 
connected to verify the calibration procedure. Before any 
vector network analyzer system utilized for measuring the 
characteristic of devices can be used, it must be 
calibrated to remove any systematic errors in the 

15 measurement hardware. Accordingly, conventional calibration 
techniques at microwave or millimeter wave frequencies 
demand the fabrication of precise short circuits, open 
circuits and matched loads which need to be fabricated to 
exacting tolerances, not easily achieved. Moreover, such 

20 methods are inhibited by dispersion in the transmission 
medium, for example, microstrip, slotline, coplanar line, 
coplanar strips, etc. in which the measurement is 
performed. These dispersion characteristics are often 
unknown in millimeter wave integrated circuit devices 

25 thereby further compounding measurement difficulty. Also, 
the integrity of the operation hinges very critically on 
whether ideal electrical characteristics are reproduced at 
the launch plane for every calibration and verification 
device and for the device under test. Thus, conventional 

30 techniques do not furnish a convenient method for assessing 
the uncertainties due to variations in such 
characteristics „ 

A typical prior art calibration technique is disclosed 
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in an article entitled, "Calibration of an Automatic 
Network Analyzer Using Transmission Lines of Unknown 
Characteristic Impedance/ Loss and Dispersion", by E. F. 
daSilva et al. appearing in The Radio and Electronics 
5 Engineer , Vol. 48, No. 5, pp.. 227-234, May, 1978. 

It is an object of this invention, therefore, to 
provide an improved method for calibrating electronic 
measurement apparatus. 

It is another object of this invention to provide an 
10 improved method for calibrating network analyzers for 
evaluating conventional and integrated circuits operating 
at microwave and millimeter wave frequencies. 

It is still another object of this invention to 
provide an improved method for calibrating network 
15 analyzers by employing measurement redundancy to enhance 
stability, convergency, and accuracy of the measurements. 

It is a further object of this invention to provide a 
method for calibrating vector network analyzers by 
utilizing relative characteristics among standards as 
20 opposed to their absolute characteristics. 

It is yet a further object of this invention to 
provide a method for calibrating vector network analyzers 
which yields a direct and convenient measure of calibration 
quality. 

25 It is still a further object of this invention to 

provide a self-verifying method for calibrating vector 
automatic network analyzers. 
Summary 

Briefly, the foregoing and other objects of this 
30 invention are provided by a method which determines the 
systematic errors associated with millimeter wave vector 
network analyzers utilized to measure complex reflection 
and transmission coefficients or impedance characteristics 
of components and devices used in microwave/millimeter wave 
35 waveguides or transmission lines. The method uses 
redundancy of measurements and relative characteristics to 
furnish information to indicate the level of non- 
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correctable corruption entering into the calibration 
procedure due to random errors. Also, as a part of the 
calibration procedure according to the present indention, 
there is provided a determination of the propagation 
5 properties of the transmission medium such as loss and 
dispersion. 

Ratios of the respective lengths of transmission line 
offsets -coupled to a measurement port of a vector network 
analyzer provide a means to generate independent 

10 measurement conditions without the need for independent 
high quality standards. The method is thus not restricted 
by dispersion, losses or the lack of a plurality of high 
quality standards. 

In the illustrated embodiment detailed below, the 

15 calibration procedure is controlled by a microprocessor and 
involves storing data for measurements taken via the 
network analyzer for a plurality of frequencies swept over 
a frequency range for six different circuit elements 
coupled to the measurement port of the vector network 

20 analyzer, namely, a load, an offset load, a doubly offset 
load, a short, an offset short, and a doubly offset short. 
Neither the offset nor the load need to be known a priori. 
They are determined by the calibration algorithm which also 
furnishes the usual error terms i.e. source match, 

25 directivity, frequency response tracking, etc. The value 
of load and/or offset as calculated by the procedure can be 
compared against an actual measurement of that component by 
the network analyzer to derive verification of the 
calibration . 

3 0 Brief Description of the Drawings 

The objects of the present invention and the attendant 
advantages thereof will become apparent by reference to the 
following drawings, wherein: 

.Figure 1 is an electrical block diagram of a vector 

35 network analyzer system utilized in practicing the method 
according to the subject invention; 

Figures 2A and 2B are electrical schematic diagrams 
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illustrative of standard error models for one-port and two- 
port devices; 

Figure 3 is a polar representation indicating the 
corruption factor for a typical calibration - run utilizing 
5 the apparatus shown in Figure 1; 

Figure 4 is a polar representation indicating the 
corrruption factor obtained when time for instrument warm- 
up was not adequate; 

Figure 5 is a polar representation indicating the 
10 corruption factor for a calibration procedure when using a 
sliding load to produce offset load conditions; 

Figure 6 is a unit circle diagram helpful in 
understanding the method according to the subject 
invention; 

15 Figure 7 is a polar representation of the magnitude 

and phase of an offset short calculated from the procedure; 

Figures 8A and 8B are characteristic curves 
illustrative of the magnitude and phase of a source match 
error coefficient obtained by the subject invention; 
20 Figures 9A and 9B are characteristic curves 

illustrative of the magnitude and phase of a reflection 
tracking coefficient obtained by the subject invention; 

Figures 10A and 10B are characteristic curves 
illustrative of the magnitude and phase of a directivity 
25 error coefficient obtained by the subject invention; 

Figures 11A and 11B are characteristic curves 
illustrative of the standing wave ratio and return loss 
representation of a reflection coefficient obtained by the 
subject invention; 
30 Figure 12 is a polar representation of the impact of 

corruption on a reflection coefficient obtained by the 
' subject invention; 

Figure 13 is an illustrative curve of the impact of 
corruption shown in a log magnitude format on a reflection 
35 coefficient obtained by the subject invention; and 

Figure 14 is a polar representation similar to Figure 
12 and being illustrative of the impact of corruption on an 
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offset short connected to the apparatus of Figure 1. 
Detailed Description 

Referring now to Figure 1, there is shown a typical 
vector network analyzer system for evaluating millimeter 
wave circuit components and including a vector network 
analyzer 10 comprising, for example, a Hewlett-Packard 
HP8510 network analyzer. The vector network analyzer 10 is 
coupled to a printer 12 and a plotter 14 by means of a 
digital bus 16 designated a system bus. The network 
analyzer 10, moreover, is controlled by a system controller 
18 comprising a microprocessor. The microprocessor 18 is 
coupled to the vector network analyzer 10 through a second 
digital bus termed a public bus which is also connected to 
a 'local oscillator source 22 and a test signal source 24 
15 whose frequency and power output is controlled by the 
microprocessor. 

The system as shown comprises a two port network to 
which is coupled a device under test (DUT) 26. One port, 
port 1, is coupled to a millimeter wave frequency signal 
through an amplifier 28 coupled to the test signal source 
24 which has been multiplied by a factor of X2 or X3 by a 
frequency multiplier 30 which is then connected to port 1 
through a signal isolator 32 and a pair of millimeter wave 
signal couplers 34 and 36. The second measurement port, 
25 designated port 2, comprises an output port which is 
coupled back through the vector network analyzer 10 through 
a third millimeter wave coupler 38, bandpass filter 40, an 
isolator 42 and a signal mixer 44 The mixer 44 is coupled 
back to the local oscillator source 22 by means of an 
30 amplifier 46. The output of the mixer 44 couples back to 
the vector network analyzer 10 by way of an amplifier 48. 
The multiplied output frequency of the test signal source 
24 is also fed to a pair of signal mixers 50 and 52 which 
are respectively coupled to the couplers 34 and 36 via the 
35 isolators 54 and 56. Mixers 50 and 52 also receive a local 
oscillator input from the local oscillator source 22 by way 
of amplifier 58. The output of the mixers 50 and 52 are 



20 
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coupled back to the vector network analyzer 10 by way of a 
pair of amplifiers 60 and 62. A typical one port 
calibration procedure involves connecting a plurality of 
standards to port 1 whereupon the system operates as a 
5 reflectometer . 

The system controller includes software which sets the 
test signal source 24, then the local oscillator source 22 
and when the system has settled makes the measurement and 
plots the data in a user-prescribed format on a CRT screen, 

10 not shown, associated with the vector network analyzer 10. 
This sequence is repeated for each of 51, 101, 201 or 401 
points over a selected frequency range in a well known 
manner for the test apparatus represented in Figure 1. 

Conventional one port calibration methods involving 

15 vector network analyzers such as shown in Figure 1 demand 
the fabrication of precise short circuits and matched 
Loads, all to exacting tolerances often beyond the reach of 
current state-of-the-art monolithic fabrication techniques. 
This also applies to known two-port calibration methods. 

20 Further, such techniques are inhibited by dispersion in the 
transmission medium, e.g., microstrip, slot line, coplanar 
line, coplanar strips, etc. These dispersion 

characteristics are often unknown in the integrated 
circuit environment, thereby further compounding 

25 measurement difficulty. Also, the integrity of a 
calibration of the vector network analyzer hinges very 
critically on whether identical electrical characteristics 
are reproduced at the launch plane (Port l)for every 
calibration and verification device and for the device 

30 under test. Thus, state-of-the-art techniques do not 
furnish a convenient method for assessing the uncertainties 
due to variations in such characteristics. Unfortunately, 
integrated circuit technology is very prone to def i'ciencies 
arising from launch plane integrity, thus leading to poor 

35 confidence in the quality of calibrations. Finally, the 
conventional approach of verifying a calibration by 
inserting a separate verification device assumes good 
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launch plane integrity for the verification standard. This 
has been found to be a poor assumption for the integrated 
circuit environment. 

The present invention overcomes these limitations by 
5 capitalizing on the strengths of integrated circuit 
fabrication. For example in such fabrication, it is easier 
to hold parameter ratios to tight tolerances rather than 
the discrete values of the parameters themselves. This 
invention, as will be shown, is well suited for such a 

10 condition because, as the following example will show, it 
is based on the use of ratioed lengths of transmission line 
offsets coupled to the measurement port, i.e. Port 1 of 
Figure 1. The algorithm to be described herein furnishes 
the transmission line's loss and dispersion characteristics 

15 as part of the calibration procedure, a decided advantage 
for integrated circuit evaluation. This invention also 
reduces the need for simultaneously maintaining extremely 
tight tolerances over a range of calibration device types. 
Instead, it can be applied based on a priori knowledge of 

20 only one standard type, e.g. a short circuit, so that only 
the assumed characteristics of this type of standard must 
be met in production. All the other standards are 
determined relative to this one standard (short), as part of 
the calibration procedure, and therefore can have an 

25 arbitrary tolerance. For these standards, however, the 
technique relies on the reproducibility of the integrated 
circuit environment. 

Lastly, this invention integrates calibration and 
verification to produce three desirable outcomes: (1) a 

30 measure of the corruption that entered into the calibration 
without the necessity of completing the process, (2) a 
quantitative display of the uncertainty limits for the 
calibration/verification procedure, and (3) a basis for 
inserting automation into the calibration process. The 

35 first outcome allows the user to abort a given calibration 
run because of the poor quality of the acquired raw data 
without ever having to compute the calibration 
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coefficients. The second gives the user an immediate 
estimate of the quality of the calibration, and the third 
provides a natural basis for artificial intelligence 
automation of the procedure so that high throughput, vital 
for reducing integrated circuit costs, is possible. 

As in other known prior art calibration techniques, 
the method of the subject invention exploits the usual 
transformation property of a terminated transmission line. 
More specifically, standard transmission line theory 
predicts that a load (L) with an inherent reflection 
coefficient of T will transform to a reflection 
coefficient T L as: 



(i) 



when measured through a length {£) of transmission line. 

15 In general, all the Ts are complex numbers, i.e. having 
both a magnitude and a phase. As defined above, z = e 22rX , 
where TT is the complex propagation factor for the 
transmission line, and T= <*+ j& , where < is the line's 
attenuation constant and is its phase constant. 

20 with the foregoing in mind, this invention is based on 

well known standard error models represented by signal flow 
graphs (Figures 2A and 2B) with coefficients expressed as 
frequency response tracking E R or E T , directivity E D , 
source match or load match E s or E L and isolation E x . 

25 These error terms are determined by calibration, and are 
alternatively called calibration coefficients or error 
coefficients. In addition to these error coefficients the 
propagation factor z of the transmission line and the 
reflection coefficient T L of at least one other standard 

30 such as a load must be determined. Thus, there are five 
unknowns and, therefore, five independent measurements are 
needed. 

The use of offset lengths of transmission line offers 
a convenient way to generate independent measurement 
35 conditions without the need for independent standard types 
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For example, adding an offset length (i ) of transmission 
line to a short circuit ( T short = -1) creates a new 
measurement of an "offset short"/ provided the larger 
offset length is less than a half wavelength at the highest 
5 frequency of interest. It -must be less than a quarter 
wavelength if the known standard is a short (open) and the 
unknown is an open (short). In the present invention 
technique, two offsets of lengths (J.) and {2 I ) are used. 
The ratio of the lengths, i.e, 2:1, is of greater 
10 importance to the technique's success than the actual value 
of U). 

If these offset lengths are combined with the short 
circuit/ three measurement conditions are generated which 
can be concisely represented by the expression: 

where C Sn is the raw or measured data value of the 
reflection coefficient T=-l/ n=0 corresponds to a short 
n=l corresponds to an offset short, and n-2 corresponds to 
a doubly-offset short. The flow diagrams of Figures 2A and 
20 2B contain the prescription for obtaining the expressions 
for Equations 2 and 3 as is well known in the art. 

A similar expression C Ln can be developed for an 
arbitrary load that is different from the short and its 
combination with the same n=l and n= 5 2 offsets, namely, 

2^- £ sTL 



Expressions (2) and (3) thus provide equations for CgQ, 

C S1' C S2' and C L0' C L1' C L2 wh i c h can be used for 
determining the five above identified unknowns of E R , E D , 
30 Eg, z, and T^. However, the six equations may be used to 
eliminate all the unknowns from the measurement set, so 
that a relationship can be developed relating to the raw or 
measured data values only. For the particular case of the 
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instant example, this relationship is: 



25 



K u - ^ (4) 



where 



K u » " C " (5) 



vc zo 



io K* = 







Cuz- 


Cu 


Cu - 


Cu> 




Csa 


Cuo - 


Cso 


Cs2 - 





(6) 



So 



(7) 



Equation (4) forms the basis for defining a corruption 
factor, K Cor as 



If the errors occurring during calibration were solely 
15 systematic, then Eq. (4) would be satisfied exactly and 
K Cor='°- 

In practice, however, non-systematic errors arising 
from noise, drift, cable stability, connector 
repeatability, spurious responses, etc., will generally 
20 cause K Cor to have a non-zero value. When a calibration 
was performed on a vector network analyzer system as shown 
in Figure 1 and including a Hewlett-Packard (HP) model 
R11643A Ka-band extender test set in conjunction with a HP 
8510A vector network analyzer, the distribution of K Cor 
over a 26.5 GHz to 40 GHz frequency range on a polar 
display as shown in Figure 3 directly indicates the extent 
to which the integrity of the calibration is compromised. 
The farther away from the origin the distribution, the 
poorer the integrity of the calibration. 
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K Cor shows very strong correlation with non- 
correctable errors. For example, Fig. 4 shows K Cor for a 
calibration performed on the same vector network analyzer 
system as shown in Figure 1 immediately after power turn- 
5 on, i.e., without allowing sufficient time for the 
instrument to stabilize. Figure 5 illustrates the effect 
of connector repeatability on the corruption factor and is 
distinguished from Fig. 3 in that the offset loads were 
produced by using a micrometer adjustment on a sliding load 

10 instead of by way of waveguide shims. 

Unlike conventional methods, information about 
corruption in the raw data is available by this invention 
immediately after data acquisition and does not require a 
calibration i .e . , determination of the error terms 

15 referred to above, to be completed. Thus, a user may 
reject an excessively corrupted data set by aborting a 
calibration run early in its cycle. In an automated 
environment, measurement throughput will thus be increased 
because of the early detection and rejection of corrupt 

20 data sets. 

At this point, redundant measurement is exploited to 
obtain an algorithm with superior numerical stability and 
convergence properties even at relatively high corruption 
levels, a distinct advantage and improvement as it relates 

25 to use in connection with millimeter wave integrated 
circuit testing. The calibration method produces stable 
outputs, e.g., offset z, load reflection T^, and the error 
coefficient terms, even though its inputs, namely the raw 
data, have fluctuations arising from corruption. Numerical 

30 stability is achieved by anticipation of the distribution 
of the raw data on a polar diagram plot and their 
susceptibility to corrupting influences from signal-to- 
noise considerations and recognizing that the calculations 
involve ratios of vector differences between the points 

35 that represent the raw data. Thus, vector differences 
between points that are closely spaced relative to the 
fluctuation in the position of those 'points will be very 
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sensitive to those fluctuations and must be avoided. Using 
redundancy is key to providing the flexibility for such 
avoidance. Convergence is manifested by the agreement of 
the determined quantities with their actual or true values. 
5 The equations represented by (2) and (3) can be used 

to eliminate all the E's and T L , so that the relationship 
between the raw data and the offset z can be expressed as: 

B * _ c 2 + i = O (9) 

The constant C, in Equation 9, is a function of the raw 
10 data C Sn and C Ln and can be expressed as such, in a number 
of mathematically equivalent forms. However, because the 
raw data is corrupt, certain functional forms are preferred 
because they are less sensitive to corrupting influences. 
Accordingly, C is expressed as: 



20 



where 



and the other constants have been defined previously. The 
expression for C in Eq. (10) has been derived in part by 
using Eq. (4) and the knowledge that K L as defined in Eq. 
(5) is more susceptible to corruption when T L is a load 
than are either Kg or K J| . This is one example of how 
redundant measurement is exploited as asserted above. 
Equation (9) can be solved to obtain the two roots, z + and 
25 z_, where 



(12) 



A number of important observations can be made regarding 
the roots. From the theory of quadratic equations, 
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(13) 



and 



2^ + 5l 



= C 



(14) 



10 



15 



20 



25 



Equation (13) suggests that the two roots must be phase 
conjugate, or, more generally., arg (z + ) + arg (z_) = 2mYT, 

where m * 0,1,2,3, Also, if the magnitude of any one 

of the roots is greater than unity, the other root must 
have a magnitude less than unity. These conditions are 
illustrated in Fig. 6. According to Eq. (1), the root that 
corresponds to attenuation, due to line loss, has a 
magnitude greater than unity, and is therefore the correct 
root. In the absence of line loss, the two roots must have 
unit magnitude and be complex conjugates of each other. 
For this case, Eq. (14) implies that C is real. Thus, the 
imaginary part of C is a sensitive indicator of 
transmission line losses and also corruption. These 
conclusions, however , are not obvious from equation 10 for 
the expression for C. 

Once z is determined, its reciprocal is the 
transformation factor for the single offset according to 
Eq. (1), and, in particular, -z"" 1 represents the reflection 
of the offset short as displayed in Fig. 7, which 
represents the magnitude and phase of the offset short 
after the correct root is chosen. 

The source match error coefficient E s is next 
calculated from the previously determined z (Eq. 13), K c 
(Eq. 8), and the formula 




(15) 



derived from Eq. (2). The result is displayed in Figs. 8a 
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and 8b which disclose the magnitude and phase of the source 
match. 

Eq. (2) may be manipulated to derive the reflection 
tracking error coefficient E R which can be expressed as: 

5 

E R = 0+ g.^f»+ ^(C+t-Cso) (16) 



15 



and is shown in Figs. 9a and 9b. 

Next, the directivity error coefficient, E D , may be 
10 obtained from Eq. (2) with n = 0, as: 

E» - C S6 + 1^ (17) 

and is illustrated in "Figs. 10a and 10b. 

Once the error coefficients E s , E R and E D are 
determined, the reflection coefficient T L of the load can 
be determined from Eq. (3) with n = 0, as: 

r c «= Ccuo-e^ — (18) 

The load reflection is shown in the standing wave ratio and 
return loss curves of Figs. 11a and lib, respectively. 

Up to this point, all of the primary parameters have 
been determined, namely, the errors terms, the reflection 
coefficient of the load T L , and the offset z, and K Cor has 
provided a measure of the corruption in the calibration 
process. However, the relationship between the displayed 
value of corruption, K Cor and confidence in the calibration 
25 procedure has not yet been established. 

Because there are six measurement conditions, the 
obvious verification strategy would be to use any five of 
them to determine the unknowns identified above, and the 
sixth relation to verify the calibration. However, such an 
algorithm is practical only for highly sanitized caw data. 



20 



30 
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Greater numerical stability of the algorithm results from 
developing the error terms described above and using the 
approach now to be described. 

Any number of the six measurement conditions may be 
5 used to assess the effect of the corruption factor, K Cor on 
the calibration. Because the reflection coefficient of the 
load T L and offset z are outputs of the procedure, the 
offset short and the load measurement are recommended for 
the assessment, thus allowing the effect of K Cor to be 

10 determined at the two extremes of the dynamic range of the 
measuring apparatus. 

Equation (4) becomes a starting point for making 
assessments of corrupting influences on all terminations. 
In particular, it can be used to solve for the raw data of 

15 the load termination, i.e., C L0 , in terms of the other raw 
data. The result is 

C u + v- __K s C So 

r L0 _ C Cl * : (19) 

| + {C.X.Z - Cu ) is 

(c L2 - c«) Ks 

The equality of Eq. (19) holds only in the absence of non- 
systematic (i.e., corrupting) influences. The right-hand 

20 side of equation 19 may therefore be viewed as an alternate 
representation of C L0 , which approaches the actual measured 
raw coefficient in the limit of zero corruption. Thus, 
both sides of Eq. (19) can be substituted in Eq. (18) to 
determine the corresponding F L . On the HP 8510A, vector 

25 analyzer 10 (Figure 1) the value of T L corresponding to 
the right-hand side of Eq. (19) is stored in memory and the 
actual load reflection (left-hand side) is viewed as data 
after the correction is applied. This procedure is 
executed with the load connected . 

30 A very sensitive indication of the impact of 

corruption is now obtained by dividing the two values of 
Pl* In the absence of corruption, the resulting display is 
a single point located at (1, 0) on a polar display. 
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Corruption will cause "smearing" of the .point as Figure 12 
illustrates for an actual load verification; Pig. 13b 
presents the data of Fig. 12 in log magnitude format, and 
Fig. 11 shows the actual return loss. The agreement is 
5 within ± 0.5 dB for a return loss in excess of 35 dB over 
the entire Ka band. 

The same procedure could be applied to an offset 
short. However, an alternate method is available based on 
the calculated offset z since the offset short can be 
10 represented by the left and right hand sides of the 
eguality: 

= Cst - £•& (20) 

( Cs, - E^E^+E* 

which holds true when corruption is absent. Alternatively, 
the value of the offset short can be evaluated by 

15 substituting its raw measurement C S1 and the computed error 
terms (Egs. 15 through 18), into Eq. (2) with n « 1, and 
comparing the resulting estimate with the previous 
representation. On dividing one by the other, a single 
point results (Fig. 13) in close proximity to (1, 0). In 

20 practice, the trace ratios are always less than ±0.005 dB 
(0 dB in a perfect system) in magnitude, and the angular 
elevation from the real axis is less than ±0.8° (0° in a 
perfect system). Better phase stability will be obtained 
if the ambient temperature is controlled. In the same 

25 manner, a comparison can be made for a doubly offset short 
by testing for the equality as represented by the 
expression: 



_ V J- / 

C Csa - Es )E* + E« 

and for a simple short with no offset it can be determined 
30 by testing the equality: 

_j - Cso ~ &3> ^ (22) 

C Cso - EtJ Er. - E-t?. 
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Thus the corrupting influences can be determined by testing 
the equalities as expressed in equations (20/ 21 and 22). 

Comparing the *0.5 dB uncertainty (Fig. 12) for the 
load measurement with the ^0.005 dB uncertainty (Fig. 13) 
5 for the short indicates that corruption has less influence 
on the short than on the load. This is to be expected if 
the chief corrupting influence is connection interface non- 
reproducibility . In practical measurement situations, 
discontinuities at the interconnect will produce small 

10 reflection amplitudes that will have a more pronounced 
influence on the low-level reflection emanating from a load 
than on the much larger reflection signal returned by the 
short. Nevertheless, this verification step for the offset 
short is strongly recommended for detecting gross operator- 

15 induced flaws during calibration. 

Although an algorithm for assessing the effect of K Cor 
on the calibration has been described, this specification 
is incomplete without a comparison of the calibration's 
outputs with the true values of the corresponding 

20 quantities, e.g., how close the offset is to its "true" 
value. In the foregoing example where a calibration is 
implemented in WR-28 rectangular waveguide using the 
dominant TE^q mode, standard waveguide theory may be used 
to predict the true value of offset at each frequency. For 

25 an actual calibration run utilizing a system such as shown 
in Figure 1, the calculated offset values are 104.1° at 
26.5 GHz and 20° at .40 GHz for the particular offset length 
(1.956 mm for the smaller offset) used, while the 
corresponding measured values (Fig. 5) are 103.89° and 

30 19.3°, respectively. In actual experience, measured 
deviations were smaller than 3° even when K Cor was 
relatively large. 

Thus, the method of the subject invention involves: 
(a) acquiring and storing raw data for the measurement 

35 conditions involving at least one known standard such as a 
short and its associated combinations of an offset short 
and doubly offset short and an unknown standard differing 



WO 89/04001 



PCT/US88/03772 



from the known standard such as a load, offset load, and a 
doubly-offset load; (b) determining, by means of software 
the corruption factor K Cor in accordance with equation (8) 
and the factors of the raw data expressed in equations (5) 
5 through (7); (c) determining the offset z from equation (9) 
taking into account the factors expressed in equations (4), 
(10) and (11); (d) finding the roots of the offset z in 
accordance with equation (12) which will provide an offset 
value for z; (e) calculating the source match, E s , from 

10 equation (15); (f) calculating the reflection tracking 
coefficient, E R , in accordance with equation (16);- (g) 
calculating the directivity, E D , in accordance with 
equation (17); (h) calculating the reflection coefficient, 
37 L , in accordance with equation (18); (i) determining the 

15 corrupting influences on all terminations particularly the 
load termination C LQ in accordance with equation (19); and 
(j) and displaying the results on the network analyzer as 
set forth above. 

There are significant benefits to having a calibration 

20 procedure in which the load, IJ. , does not need to be known 
a priori. For example, one may use the device on which 
measurement is ultimately desired (and for which the 
instrument is being calibrated in the first place) as the 
"load". The procedure will then determine its reflection 

25 coefficient. Another significant example is one where the 
return port, i.e. port 2 of Figure 1 is used in place of 
"T L . In this case, the load match error term, E L , is 
determined by the method set forth above. The other two- 
port error terms may be obtained by combining this 

30 technique with other conventional methods. 

Having shown and described what is at present 
considered to be the preferred method calibration, it 
should be known that the same has been made by way of 
illustration and not limitation. Accordingly, all 

35 modifications, alterations and variations coming within the 
spirit and scope of the invention as set forth in the 
appended claims are meant to be included. 
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Claims; 

1. A method for calibrating measurement apparatus 
such as a vector network analyzer for making evaluations of 
circuits operating at microwave and millimeter wave 
frequencies comprising the steps of: 

(a) acquiring and storing data for measurement 
conditions involving an arbitrary known standard 
termination alone and in combination with a ratioed length 
pair of offsets; 

(b) acquiring and storing data for measurement 
conditions involving an unknown standard termination 
different from the known standard termination alone and in 
combination with a pair of offsets having the same length 
ratio as set forth in step (a); 

(c) determining a figure of merit of the quality of 
15 the acquired raw data; 

(d) determining the value of the offsets used in 
steps (a) and (b); 

(e) determining a plurality of predetermined error 
terms or coefficients; 

20 ( f ) determining the reflection coefficient of the 

unknown standard; and 

(g) determining the corrupting influences on said 
terminations by comparison of their calculated values as 
determined in steps (d) and (f) with their measured values. 

25 2 - The method as defined by claim 1, wherein the 

known standard comprises a short circuit, an open circuit, 
or a termination with known but arbitrary response. 

3. The method as defined by claim 1 wherein the 
unknown standard comprises a load, a short circuit, an open 
30 circuit, or other termination with unknown and arbitrary 
response different from the response of said known 
standard. 
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4. The method as defined by claim 1 wherein said 
offsets have a length ratio of 2:1. 

5. The method as defined by claim 1 wherein said 
figure of merit of step (c) comprises the corruption factor 

5 K Cor* 

6. The method as defined by claim 1 wherein said 
error coefficients include frequency response tracking, E R , 
directivity, E D , and source match, E s , and/or load match, 
E L . 

10 7 » Ti *e method as defined by claim 6 wherein said 

error coefficients additionally include isolation, E x , and 
transmission tracking, E T . 

8. The method as defined by claim 1 wherein said 
figure of merit of step (c) comprises the corruption factor 
15 K Cor whicn is expressed as: 



20 



where 



K 







ic Zo 






Cu, 




Cu - 


Cuo 




.Ccz - 


Csz 




Cut — 






Cs*. - 


Csi 



and 



Csi - C 



so 



where C SQ , C sl , C s2 , C L0 , C L1 and C L2 comprises measured 
data values of the reflection coefficients for a short, 
offset short, doubly offset short, a load, offset load, and 
doubly offset load, respectively. 

9. The method as defined by claim 8 where said step 
(d) comprises determining the offset z from the roots of 
the expression: . 

Z - C 2 +1-0 
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where 



and 

K 



10 a Cu - Cst 



10. The method as defined by claim 9 wherein said 
roots are expressed as: 



C*-4 



5 11. The method as defined by claim 10 and wherein one 

of the error coefficients of step (e) comprises the source 
match. Eg, which is determined by the expression: 

F - - » C Lz 



12. The method as defined by claim 11 wherein one of 
the error coefficients of step (e) comprises the reflection 
10 tracking error coefficient, E R , which is determined by the 
expression: 



13. The method as defined by claim 12 wherein one of 
the error coefficients of step (e) comprises the 
directivity error coefficient, E D , which is determined by 
15 the expression: 



1+ E s 



14. The method as defined by claim 13 wherein said 
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step (f) of determining the reflection coefficient T L , in 
accordance with expression: 

P L = (Co - E P ^ 



15. The method as defined by claim 14 wherein said 
step (g) of determining corrupting influences is determined 
from the expression: 

t-U T 7— ~ N^C- 

/* _ {Ci.Z - Csz) 
v-i_o — — - — L 



»s so 

•uo 

. L2 - Cui^ 



(Cu - Cm) ^ 

by using both sides of this equation to obtain estimates of 
II as set forth in step (f), and testing for the equality 
of the two estimates. 

10 16 ' The method as defined by claim 8 wherein said 

step (g) of determining corrupting influences is determined 
by testing the equality: 



17. The method as defined by claim 16 wherein said 
step (g) of determining corrupting influences is determined 
15 by testing the equality: 



" 2 Cs 2 - 



-2 = 



18. The method as defined by claim 16 wherein said 
step (g) of determining corrupting influences is determined 
by testing the equality: 

C So — Et> 



( do- e^e^ + E 
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